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Abstract — We consider the design of coding schemes for the 
wireless two-way relaying channel when there is no channel 
state information at the transmitter. In the spirit of the compute 
and forward paradigm, we present a multilevel coding scheme 
that permits the recovery of a class of functions at the relay. 
We define such a class of functions and derive rates that are 
universally achievable over a set of channel gains when this class 
of functions is used at the relay. We develop our framework 
with general modulation formats in mind, but numerical results 
are presented for the case where each node transmits using 
the QPSK constellation. Numerical results with QPSK show 
that substantially higher rates are achievable with our proposed 
approach than those achievable by always using a fixed function 
or adapting the function at the relay but coding over GF(4). 

Index Terms — Network coding, multilevel coding, two-way 
relaying, compute-and-forward 

I. Introduction 

Physical layer network coding (PLNC) or Compute and 
Forward is a new paradigm in wireless networks where each 
relay in a network decodes a function of the transmitted 
messages and broadcasts the value of this function to the 
other nodes in the network. This has been shown to provide 
significant increase in achievable rates for some networking 
problems fH, Q, Q. An example of such a problem where 
compute and forward has been shown to be effective is the 
two-way relaying problem shown in Fig. [T] Node A has data 
to send to node B and vice versa. The relay R is included to 
assist in this communication, and it is assumed that there is 
no direct link between nodes A and B. Near optimal coding 
schemes have been designed for this problem for the case 
where there is no fading in the channel in 0, Q, (4) and for 
the case when there is fading but the transmitter has perfect 
channel state information in 0. 

In this paper, the complex channel coefficients Ha and Kb 
are assumed to be perfectly estimated at each receiver but 
unknown to each transmitter. For this scenario, the authors in 
(6) introduce a scheme called denoise-and-forward which uses 
channel dependent denoising functions at the relay to minimize 
the symbol error probability. The relay chooses denoising 
functions so that the distance profile for constellation points 
with different labels is optimized. This improves the symbol 
error rate for transmissions between nodes A and B, however, 
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Fig. 1. Diagram of two-way relaying channel with PLNC. 



denoising is performed purely at the symbol level. There is no 
natural extension to include error correction at the relay. 

Recently, a scheme called compute-and-forward, which al- 
lows both adaptation of decoding functions and error correc- 
tion at the relay has been presented in Q. In this scheme, 
the relay decodes an integer combination of the transmitted 
codewords, where the integer combination is adapted accord- 
ing to the channel gains. They show that such a scheme can be 
implemented using nested lattice codes to take advantage of 
the duality between modulus arithmetic on prime order fields 
and the modular operations of lattice decoding. Their scheme 
requires the construction of infinite dimensional lattices which 
is not practical. The results in [7 ) are extended in a remarkable 
way in [ 8 ], where an algebraic framework is provided to design 
lattices over principal ideal domains. However, their proposed 
coding scheme is also based on large dimensional lattice codes. 

In this paper, we propose a compute and forward scheme 
based on multilevel coding (MLC). Unlike the coding schemes 
in 0, 0, our proposed scheme does not result in a lattice 
code and uses only linear codes over prime fields (for example, 
binary linear codes) and, is hence, practical. Yet, it facilitates 
error correction for a larger class of decoding functions than 
those proposed in (7). This is because the class of functions for 
our scheme is derived from the large set of non- singular square 
matrices over ¥ p in place of the set of non-zero elements in 
large prime order fields. To the best of our knowledge, such 
an idea of using multilevel coding and exploiting the linearity 
over the prime field to adaptively decode linear functions of 
transmitted codewords is new. Another important contribution 
in this paper is that our proof for the achievability of rates with 



the proposed multilevel coding scheme requires a non-trivial 
extension of the proof of achievability of rates for multilevel 
coding for the point to point case. 

This paper is organized as follows. The key elements of 
the problem are outlined in Section [TT| Our proposed solution 
is detailed in Section [Hi] The rate which can be achieved by 



our scheme for the MA stage is given in Section IV These 
rates are numerically determined for an example where nodes 
A and B transmit using a QPSK constellation in Section |V| 



Key results are reiterated in Section VI 



II. Problem Description 

Each node in the relay network is assumed to be half-duplex, 
so communication is split into two stages, a multiple access 
(MA) stage and a broadcast (BC) stage. We assume perfect 
synchronization between the transmitters and mainly focus on 
the MA stage in this paper. 

A. Multiple Access Stage 

Nodes A and B each encode their binary messages u A 
and u B into codewords v A G C A and v B G C B where 
Ca and Cb are the codebooks used at the nodes A and B 
respectively. These codewords are mapped to sequences of 
symbols s A , s B G Q N with \Q\ — 2 £ . The relay receives noisy 
observations of the sum of these symbol sequences according 
to 



h A s A + h B s B + w R 



(1) 



where h A and h B are complex fading coefficients, and w R is 
complex additive white Gaussian noise (AWGN) with variance 

B. Adaptive Decoding at the Relay 

The main idea proposed in this paper is the construction 
of a coding scheme such that the relay can reliably decode 
some function of v A and v B for a desired set of channel 
conditions H C C 2 . Specifically, we jointly design codes 
Ca and Cb and a set of decoding functions T such that, 
for any (Ka^b) £ H, there exists / G T such that the 
relay can reliably decode f(v A ,v B ) from y . We require 
that, given the output of f{y_ A ,v B ), node A (B) must be able 
to unambiguously decode v B (v A ) using its knowledge of 
v_ A (V-b)- For a gi yen / G J 7 , we will define an induced 
codebook at the relay as the codebook corresponding to / 
i.e. 

C f ,R = {f(v A ,v B )\v A G C A ,v B e C B } (2) 

It is important to understand the structure of C/ ?j r since the 
probability of error in decoding f(v A ,v B ) from y depends 
on h A , h B , and C/ }J r. The main advantage of our proposed 
scheme is that it guarantees that choosing one codebook 
Ca and Cb at the transmitter can result in a good induced 
codebook C/ }j r for a class of functions T. More specifically, it 
guarantees C/ ?j r is a member of the ensemble of random coset 
codes which is an optimal ensemble for achieving the uniform 
input information rate for the equivalent channel between 
f(v A ^ B ) and y R for all / G T. 




Fig. 2. Block Diagram of MLC Coset Encoders for MA Stage. 



The broadcast stage is fairly standard and is identical to that 
considered in 0, (3). 

III. Proposed Scheme 
A. Multilevel Encoder 

The system model for the multilevel encoder for nodes A 
and B and the channel model for the MA stage is shown in Fig. 
12 The encoder at nodes A and B uses MLC with a different 
coset of the same linear code C used at each bit level. For a 
detailed description of MLC and achievable rates for the point 
to point channel, see 0. 

The encoder is described as it pertains to node A to 
simplify notation. First, the message u A is split into sub- 



vectors u 



l 
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which form rows of an £ x K matrix 



L —A 



(3) 



Each u A , {1, ...,£} is encoded with a linear code C with gener- 

These codewords 



ator matrix G to get codewords 7^ , . . . , 7 
from the rows of an t x N matrix T A 



A' 



U^G. Finally, a 

A ia ciuucu lu cav^il 7^* 

can be thought of as coset leaders of a random coset of 
the original linear code. We obtain a codeword of a random 
coset given by v A = © 7^, k G {1, ..,£}. The random 
coset leaders form an t x N matrix A^. The resulting coset 
codewords v A form the rows of a binary £ x N matrix 
given by 



random binary vector A \ is added to each 7^. Since each 



[x A [ll...,x A [N]}. (4) 



Thus each code C A , k G {1, will be a different coset 
of C. The kth row v A of X^ is then a codeword of C A , and 
the n t h column x A [n] is the binary address vector of the n t h 
symbol of s A . The n t h binary address vector x A [n] G F2 
maps to a symbol s A [n] G Q through the use of a symbol 
mapping function M. : ¥2 —> Q. It should be mentioned here 
that much of the intuition about the main result in the paper 
is best obtained by ignoring the fact that cosets are used at 
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Fig. 3. Effective constellation at relay for different values of ( 



each layer and simply considering the use of identical linear 
codes at each level in the MLC scheme. The coset matrix 
is included to symmetrize the effective channel at the relay 
(i.e. A a is necessary for the proofs to be correct). 

B. Adaptive Decoding at the Relay 

As mentioned previously, the goal of the proposed scheme 
is to allow the relay to decode a function of the transmitted 
codewords. If nodes A and B encode their messages as 
described, the set of decoding functions T which the relay 
can use for decoding is defined as follows. 

Define V as the set of i x t binary matrices which are 
invertible using operations over F 2 . The set of functions we 
consider is given by 



T = {f : Wi x Fi Wi 
I f(x A ,x B ) = [D A D B ] 



, D A ,D,GD}. (5) 



Therefore a given / G T is defined by some D^,D# G V 
from which the relay should attempt to decode a matrix X f r 

_ X ^ _ 

Due to the linearity of [Da,Db] and G, we can express 
the desired matrix X/ ?j r as 



given by X fjR = [D A D B ] 



= [D A D B ] 



A f,R- 



UaG 



)A A 
)A B 



(6) 



Here, we see that the matrix Xj ^ can be written in terms of 
an effective message Uj,r and coset matrix A/ }j r which can 
be computed separately based on /. Thus the rows of X/ ?j r 
are codewords from a different coset code of C. Note that / 
is applied elementwise to the sequences s A and s B . 



To illustrate the importance of choosing the decoding func- 
tion / depending on (hAjhs), consider an example with 
Q = = {M(00),M(01),M(U),M(10)} 

(i.e. QPSK with Gray Labeling). Further, let h A = e J ° A , 
h B = e j6>B , and = 6 a — b be the phase difference. 

\r^^- /T\o^f t»2 rx\ 



Consider the decoding functions fi (x A ,x B ) 



x B \, f2(x A ,x B ) = [x\ © x\,x\ © x B \. The resulting 
constellation Qr at the relay is shown for different values 
of in Fig. [3] Note that the complex coordinates of the 
constellation points are exactly the same, but their labels are 
different based on and / G /b}. When ~ 0, fi appears 
to have better performance than in terms of the distances 
between points with unequal labels. The situation is reversed 
when 6 ~ ^. This shows that the performance for a fixed 
decoding function can vary widely with 6 even when both 
\h A \ and \h B \ are large. 

In order for nodes A and B to be able to unambiguously 
decode their desired messages, the authors in [6] show that / 
must satisfy 

f{x A ,x B ) ^ f(x' A ,x B ) V^/ x' A and x B 
f(x A ,x B ) ^ f(x A ,x B ) V x B ^ x B and x A . (7) 

We call functions that satisfy this property unambiguous. 
Lemma 1: For any D^,D# G V, a decoding function 



f(x_AjX L 



is unambiguous. 



Proof: The full proof is^omitted for space, but follows 
from the invertibility of and T> B . ■ 

IV. Achievable Information Rates 
A. Achievable Rate for a Given Function 

For a given / and fixed channel gains h A and h B the 
achievable rate region is given by the following theorem. This 
theorem is the key contribution of this paper. 

Theorem 1: Choose some fixed D^,D# G V and define 



l -f,R 



f(x 



(8) 



Choose a subset S C {1, ...,£} and define S = {1, ...,£} \ S. 
Divide S into p non-empty disjoint subsets Si, ...,S p so that 
Ui=i $i = S> Let Z{, i G {l,...,p} define p i.i.d. Bernoulli 
random variables with parameter \. At last, let each row of 
X^ and X^ be encoded using a different coset of the same 
linear code C. Then there exists a linear code C of rate 7Z 
for which the relay can reliably decode X/ ?j r as long as 1Z 
satisfies 

l I{Y R -{xl R \keS}\{xl R \keS], 



1Z < _min 

<s,<s,«Si,...,«s p V 



{X k 



f,R 



Zi\keSi}Vie{l,...,p}). (9) 

2, the set of bounds described 



For the special case when t 
by ^ are equivalent to 



U < min{-I(Y R ;X R ,X R )J(Y R ;X R \X R ),I(Y R ;X R \X R ), 



I(Y R ;X R ,X R \X R i 



' Z 1, X R ' 



<Zi)}- 



(10) 



Note that 

i(Y R] x^x% ixiezi.^ezo = i^x^xl^xl 

That is, {X^ © Z x , X|, © Z x } and {X X R © X 2 R } carry the same 
information about X R and X R . 

Proof: The detailed proof is not included here because 
of space limitations. However, the key steps in the proof are 
outlined below. 

Our proof uses the standard approach of deriving upper 
bounds on the probability of error for a joint typicality decoder 
averaged over a carefully chosen ensemble of codes. The 
ensemble considered here is the ensemble obtained by using 
random cosets of the same linear code for the different sig- 
nalling levels in the multilevel coding scheme. The linear code 
is chosen from the ensemble of linear codes with randomly 
chosen entries in the generator matrix. The use of the same 
linear code in each level is an important ingredient in our 
proposed scheme since we allow the relay to freely take 
linear combinations of codewords from different signalling 
levels. However, this is also what complicates the proof. 
The ensemble used here is different from the often used 
ensemble of random coset codes used at each level in the 
multilevel coding scheme since the latter ensemble allows for 
independently chosen codes at each level. While the latter 
ensemble has been used widely to obtain achievable rates for 
MLC for the point to point channel and the multiple access 
channel, the former ensemble has not been analyzed in detail 
in the literature. The key contribution of our proof in the 
Appendix is to derive the achievable rates with the former 
ensemble with identical linear codes at each level. 

This can be accomplished since the use of the same linear 
code at each level ensures that for each / G J 7 , Cf R , k G 
{1, £} is a member of the ensemble used at the transmitters. 
The main complication that arises from this is that the pairwise 
independence assertion that is required in typical channel 
coding proofs |[T0l does not hold for certain classes of error 
events. Particularly, it is possible for the relay to correctly 
decode some rows of Xj^ while others may be in error. We 
handle this by splitting the union bound for error probability 
into separate classes of error events which are conditionally 
pairwise independent. 

The bound for the £ = 2 case can be derived by letting 
<S,<S,<Si,<S2 ^ take the following values respectively. 

{5 = {1,2},5 = 0,5i = {1},5 2 = {2}} 

{5 = {1},5 = {2},5 1 = {1}} 

{<S = {2},<S = {1},<S 1 = {2}} 

{S = {1,2}, 5 = 0, Si = {1,2}}. (11) 

Notice that the first three terms in ( fT0| ) are also required by 
the proof for the general MAC channel. The last bound is a 
result of the requirement that each signaling level uses a coset 
of the same linear code. ■ 

B. Universally Achievable Rate 

We say that a rate 1Z is universally achievable over the set 
% C C 2 if there exists a fixed linear code C of rate 7Z and 



coset matrices A a and A# such that for every (h A , Kb) G 
the relay can reliably decode X/ 5 # for some / G T. That is 
some X/ ?j r can be decoded with arbitrarily small probability 
of error in the usual information-theoretic sense. The main 
result in this section is the following theorem. 

Theorem 2: For a fixed / G T and (fiA^hs), define 
Tlf(hAjhB) as the supremum of rates satisfying ^ where 
£f,R ~ K^a^b)- For any finite set of channel gains, 
HcC 2 , any rate 1Z such that 

1Z < min max 7Zf(hA,hB) (12) 

(h Aj h B )en fer J 

is universally achievable. 

Proof: The above result follows from the fact that the 
MLC scheme with transmission rate 1Z is universal in the 
sense that the induced codebooks are simultaneously optimal 
for decoding any function / G T if 1Zf(hA, h B ) > 7Z- 
This permits the relay to decode any function / for which 
1Z < 7Zf(hA,hB) since Hf(hA,h B ) corresponds to the 
information rate corresponding to the uniform i.i.d distribution. 
Hence, for every {Ha-, h B ) G H we can choose a function / at 
the receiver such that 7£/(/ia, h B ) > TZ, and the result follows. 
The proof for the existence of a single coset code {C, A^, A B } 
which is good for the finite set H is omitted for space. ■ 
Note that in order for this problem to be practically inter- 
esting, the set H should be meaningfully defined. 

V. Numerical Results 

A. Numerical Results for QPSK 

As an example, consider the case where nodes A and 
B transmit symbols from a QPSK constellation with Gray 
Labeling. Fig. |4(a)| shows a plot of the achievable information 
rate £7Zf(hA,h B ) as given in flQ| for each function / G T 
dependent on the phase difference = 9a — Ob for an SNR 
of 7 dB. H is the set of channel gains 

H = {(h A , h B )\h A = e? eA ,h B = d 9 *} (13) 

where Qa,®b £ {0, ^, 2ty} for a finite integer m. Thus 
\H\ is finite but approximates a the selection of any value of 
9a and 9b arbitrarily closely. The dotted line indicates the 
universally achievable rate in bits per complex symbol for the 
proposed scheme which satisfies Theorem 2 for T~L. Note that 
different functions provide the best performance for different 
values of 9 which reiterates the substantial benefit of decoding 
adaptively. Notice that a small increase in rate makes reliable 
decoding impossible for any / G T for a significant range of 
9\ however, there are many (h A ,h B ) T~L such that 3f G T 
for which reliable decoding is possible. 

B. Coding over GF(4) 

It is interesting to use this QPSK example to compare our 
MLC scheme to the case where nodes A and B encode their 
a data in F4 using a linear code Cqfa of rate IZgfa- The 
relay uses the set of decoding functions Tqf^ corresponding 
to linear combinations of codewords in F4 of the form 

v R = f(v A ,v B ) = av A ®Pv B , a,/3GF 4 \{0}. (14) 
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Fig. 4. Numerical Results 



The value of J(Yr; /(Xa^Xb)) for each possible / G J 7 gf4 
is plotted as a function of in Fig. |4(b) with an SNR of 7 dB. 
Again the dotted line represents universally achievable rate for 
the U in {B). 

C. Comparison of Proposed Techniques 

These numerical results illustrate that the proposed MLC 
scheme facilitates better decoding flexibility at the relay than 
coding over F 4 for this example. In fact, an analysis of 
these functions based on the labeling of points in Qr shows 
that J~gfa C T. However, this improved flexibility comes 
at the cost of additional rate constraints on each / G T. 
The thick dashed line in Figs. |4(a)| and 4(b) represents the 
rate which is achievable if the relay decodes using some / 
which is equivalent to the componentwise xor operation for 
multilevel coding or finite field addition for F4. The difference 
between these curves illustrates the effects of the additional 
rate constraints imposed by ([9]). In Fig. |4(a)| the last term 
KYr-X^X^X^ © Zx,X\ Z{) in (H^is dominant if 
6 « 77 for determining the achievable rate for this function. In 



Fig. |4(b)| we see that this term does not need to be satisfied if 
nodes A and B use a linear code in F4. 

The universally achievable rate for the H in ( [T3] ) (i.e. the 
constant value given by the dotted line in Figs. 4(a) and 4(b) ) 
is plotted as a function of SNR in Fig. |4(c) for the cases 



where the relay uses T or Fgfa- This value asymptotically 
approaches 1.5 bits per symbol for coding over F 4 . From Fig. 
|4(b)| this appears to occur because Fqfa does not provide the 
relay with a decoding function which works well when 6 « | . 
This represents an extreme case, because the event \Ha\ = 
\Hb\ occurs with probability zero for many random fading 
processes. However, this illustrates that for PLNC it is possible 
for the universally achievable rate to be limited by specific 
(h>A,h>B) £ T~i even if each \Ha\ and \hs\ is large. While the 
proposed scheme is better than choosing only one decoding 
function, it still does not entirely eliminate the interference that 
results from the signal constellations not aligning perfectly at 
the receiver. Hence, the achievable rate in the high SNR regime 
will be limited by this interference in addition to the limitation 
imposed by a finite-sized constellation. 



VI. Concluding Remarks 

In this paper, we have proposed a coding scheme based on 
MLC for compute and forward or PLNC for the case when the 
channel is perfectly estimated at each receiver but unknown 
to each transmitter. We showed that MLC allows for decoding 
of a set of functions of the transmitted messages and the 
relay can choose one function from this set depending on the 
channel coefficients. In Theorem 1, we obtained an achievable 
rate for a fixed decoding function and channel realizations. In 
Theorem 2, we obtained a numerically computable expression 
for the universally achievable information rate over a set of 
channel realizations. Numerical results for QPSK suggest that 
the proposed scheme significantly outperforms the use of a 
fixed decoding function with binary linear codes and is better 
than using linear codes over F 4 . 
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